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O 

o 

(N 



A. A. Schekochihm/'E S. C. Cowley, 1 ' 2 S. F. Taylor, 1 G. W. Hammett, 3 J. L. Maron, 4 - 5 and J. C. McWilliams 6 

'Plasma Physics Group, Imperial College, Blackett Laboratory, Prince Consort Road, London SW7 2BW, UK 
2 Department of Physics and Astronomy, UCLA, Los Angeles, CA 90095-1547 
s Plasma Physics Laboratory, Princeton University, Princeton, NJ 0854-3 
4 Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 
5 Department of Physics and Astronomy, University of Iowa, Iowa City, IA 52242 
Department of Atmospheric Sciences, UCLA, Los Angeles, CA 90095-1565 
(Dated: February 2, 2008) 

We consider the problem of incompressible, forced, nonhelical, homogeneous and isotropic MHD 
turbulence with no mean magnetic field and large magnetic Prandtl number. This type of MHD 
turbulence is the end state of the turbulent dynamo, which generates folded fields with small- 
scale direction reversals. We propose a model in which saturation is achieved as a result of the 
velocity statistics becoming anisotropic with respect to the local direction of the magnetic folds. The 
model combines the effects of weakened stretching and quasi-two-dimensional mixing and produces 
magnetic-energy spectra in remarkable agreement with numerical results at least in the case of a 
one-scale flow. We conjecture that the statistics seen in numerical simulations could be explained 
as a superposition of these folded fields and Alfven-like waves that propagate along the folds. 
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In this Letter, we consider what is perhaps the old- 
est formulation of the MHD turbulence problem dating 
back to Batchelor's work in 1950 Q: incompressible, ran- 
domly forced, nonhelical, homogeneous, isotropic MHD 
turbulence described by 



<9 t u + u • Vu = vAu - Vp + B • VB 
<9 t B + u • VB = B • Vu + ijAB. 



(1) 
(2) 



The pressure p (determined from V • u = 0) and the mag- 
netic field B are rescaled by p and (47T/?) 1 / 2 , respectively 
(/? is density). Turbulence is excited by the random exter- 
nal forcing f . No mean field is imposed. We are primar- 
ily interested in the case of the large magnetic Prandtl 
number Pr m = vjr\ which is appropriate for the warm 
interstellar medium and cluster plasmas 0. Numeri- 
cal evidence suggests that the popular choice Pr m = 1 
is in many respects similar to the large-Pr m regime 0. 
Pr m 3> 1 implies that the resistive scale £ n ~ Pr" 1 ' 2 ^ 
is much smaller than the viscous scale l v . Thus, the 
problem has two scale ranges: the hydrodynamic (Kol- 
mogorov) inertial range 1$ ^> £ 3> l v ~ Re _3 ^ 4 ^o (^o is 
the forcing scale) and the subviscous range l v » I » t v . 

For a moment, let us consider the traditional view of 
fully developed incompressible MHD turbulence in the 
presence of a strong, externally imposed mean field. This 
view is based on the idea of Iroshnikov |4j and Kraich- 
nan || that it is a turbulence of strongly interacting 
Alfvcn-wave packets. This phenomenology, modified by 
Goldreich and Sridhar |(| to account for the anisotropy 
induced by the mean field, predicts steady-state spectra 
for magnetic and kinetic energies that are identical in 
the inertial range and have Kolmogorov fc~ 5 / 3 scaling. 
An essential feature of this description is that it implies 
scale-by-scale equipartition between magnetic and veloc- 



ity fields: indeed, <5uk = <5Bk in an Alfven wave. Simula- 
tions appear to confirm Alfvenic equipartition provided 
there is an imposed strong mean field Bq ;§> w rm s Q- 

In the case of zero mean field, it has been widely as- 
sumed that essentially the same description applies, ex- 
cept it is the large-scale magnetic fluctuations that play 
the role of effective mean field along which smaller-scale 
Alfven waves can propagate. However, numerical simu- 
lations of isotropic MHD turbulence do not show scale- 
by-scale equipartition between kinetic and magnetic en- 
ergies. There is a definite and very significant excess of 
magnetic energy at small scales. This is true both for 
Pr m > 1 and Pr m = 1 (Fig. P). This result persists at 
the highest currently available resolution (1024 3 , see ||). 

Let us consider the genesis of the magnetic field in 
isotropic MHD turbulence. As there is no mean field, all 
magnetic fields are fluctuations generated by the small- 
scale dynamo. This type of dynamo is a fundamen- 
tal mechanism that amplifies magnetic energy in chaotic 
3D flows with sufficiently large magnetic Reynolds num- 
bers and Pr m > 1. The amplification is due to random 
stretching of the magnetic-field lines by the velocity field. 
During the kinematic (weak-field) stage of the dynamo, 
the magnetic energy grows exponentially in time, its 
spectrum is peaked at the resistive scale, k v ~ Pr;^ 2 fc„, 
and grows self-similarly 0, 0, ■ The growth rate is of 
the order of the turnover rate of the fastest eddies, which, 
in Kolmogorov turbulence, are the viscous-scale ones. 

Although the bulk of the magnetic energy is at the re- 
sistive scale, the dynamo-generated fields are not at all 
randomly tangled, but rather organized in folds within 
which the field remains straight up to the scale of the flow 
and reverses direction at the resistive scale 0, El E^, El ■ 
One immediate implication of the folded field structure 
is the criterion for the onset of nonlinearity. For in- 
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FIG. 1: Energy spectra in simulations with Pr m = 1, 
Re a ~ 155 and with Pr m = 10, Re A ~ 45 (bold lines). The 
thin lines are our model-predicted spectra of the folded field 
component (normalized to have the same energy as the nu- 
merical spectra). 



compressible MHD, back reaction is controlled by the 
Lorentz tension force B • VB ~ k\\B 2 . This quantity de- 
pends on the parallel gradient of the field and does not 
know about direction reversals (fen ~ k v 13]). Balanc- 
ing B ■ VB ~ u ■ Vu, we find that the back reaction 
is important when the magnetic energy becomes compa- 
rable to the energy of the viscous-scale eddies. Clearly, 
some form of nonlinear suppression of stretching motions 
at the viscous scale must then occur. However, the eddies 
at larger scales are still more energetic than the magnetic 
field and continue to stretch it at their (slower) turnover 
rate. When the field energy reaches the energy of these 
eddies, they are also suppressed and it is the turn of yet 
larger and slower eddies to exert dominant stretching. 
The folded structure is preserved with folds elongating 
to the size £ s of the dominant stretching eddy. The key 
question is whether £ s can increase all the way to the 
outer scale or stabilizes just above the viscous scale [lij . 

The nonlinear suppression of stretching motions does 
not mean complete elimination of all turbulence: only 
the bb : Vu component of the velocity-gradient tensor 
leads to work being done against the Lorentz force and, 
therefore, must be suppressed. It is then natural to ex- 
pect a local anisotropization of the velocity field. In this 
Letter, we demonstrate how a simple model accounting 
for this nonlinearly induced local anisotropy can produce 
solutions that are in remarkably good agreement with 
numerically observed magnetic-energy spectra. 

The idea is to use the standard Kazantsev 91 model ve- 



locity, Gaussian and white in time, {u l (t, ~x)u 3 (f , x')) = 
S(t — t')^ 13 (x — x'), but let k 13 depend on the local direc- 
tion of the magnetic field, b = B/£?. In the Lagrangian 
frame (with local rotation transformed out), b orients it- 
self along the stretching Lyapunov direction of the flow, 
which stabilizes exponentially in time |l5j . Therefore, in 
this frame, b can be assumed to vary slowly with time. 
In the presence of one preferred direction defined by b l b J , 
the velocity correlator in k space has the following form 



K ij '(k) 



K«(fc,M)(*«- 



ki kj J 



K^\k,\^\){¥V 



fikib 3 



(3) 



where k = k/fc, fi = k • b. Let us ignore the spatial 
dependence of all quantities that vary at the flow scale 
and slower. The velocity will only enter via its gradient 
Uj = dju'\ which is now a function of time only with 
statistics (tt^ = S(t - t') J d 3 kk n k m K i3 (k). We 
can assume that b also depends on time only, because 
it will always enter via the tensor W , which varies at 
the scale of the flow (because of the folded structure of 
the magnetic field, the field's curvature is very small 3, 
H^ |. so the fast spatial variation of b is limited to sign 
reversals and cancels in b b 3 ). With these assumptions, 
the solution to Eq. (j2J) can be written as (cf. |12lllfi| ) 



B(i,x) = b(<) / d 3 k B(t,\L )e ixi ( t ' Wo \ 



where k(0, ko) = ko and 



d t B = b l b m u l m B -rjk^B, 
d t ti = b m ut 



b l b m ulb l 



(4) 



(5) 
(0) 
(7) 



Equations JSH7J) are a modification of the so-called zero- 
dimensional model of the dynamo ^t|- A closed equa- 
tion can be obtained for the joint PDF of B, k, and 
b, V(B,k,b) = <5(|b| 2 - l)<5(b • k)(4 7 r 2 fc)- 1 P(i?,fc), 
via an averaging procedure analogous to, e.g., the one 
in Ref. The magnetic-energy spectrum M(k) — 

(1/2) J °° dB B 2 P(B 7 k) is then found to satisfy 



d t M 



1 d_ 

8 1± dk 



(l + 2<7||)£r 



~dk 



(1 + 4ct_l + 10o-||)fcM 
2(cr_L + <T||)7_lM - 2r]k 2 M, 



(8) 



(1/V) / d 3 kk 2 K± , a 

= ft/*, K1 



where j± = J ' d 3 k k^K^, a± - 
(l/ 7± )/d 3 fcfc 2 K ||, k± = fc(l 
(1/2) ((^' -Vb 3 )^, «|| = (1/2)6*6^ K ij , and K ij is defined 
in Eq. In the isotropic case, = = 0, 

which gives <r± = 2/3, mi = 1/6. Equation (jHJ then 
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reduces to the standard equation for the magnetic-energy 
spectrum in the kinematic dynamo 0, Hoj . With a zero- 
flux boundary condition imposed at low k [l4|, Eq. (JSJ) 
has an eigenfunction (in the limit r\ — » +0) 



M(k) ~ k s e^K (k/k n ), 



(9) 



where TQ) is the 
1 1/2 



Macdonald function. fc r 



[(1 



2tr||)7j./16»7] , s = 2(ct ± + 2ct,|)/(1 + 2<7||), and 7 = 
( 7 x/8) [l6(ox. + an) - (1 + 2a ± + 6 CT ||) 2 /(1 + 2a,,)] . As 
magnetic back-reaction makes velocity more anisotropic, 
the values of <tj_, o~\\ drop compared to the isotropic case, 
and so does the growth rate 7 — until the dynamo is shut 
down (for a purely two-dimensional velocity, a±_ = an = 
and 7 = — 7_l/8). Thus, saturation can be achieved purely 
by anisotropizing the statistics of the velocity field. 

How do we make connection from a theory based on the 
(^-correlated model velocity to the real turbulence, which 
has a finite correlation time? The simplest prescription 
is to get finite expressions for equal-time velocity correla- 
tors by replacing the 5 function by l/r c : (u l (k)tz J *(k)) = 



J« (k) 



'(k). We take the correlation time t c 



of a given type of motions to be their "turnover time" : 
defining 7j_ and 7n analogously to k± and Km, we write 
k 2 ,Ki = C\ iTT^/i , k?,K± = C\\ ± (cr±j±y 1 I±, and 



C\\±, C|| 



jd 3 kkii ± 



1 are ad- 

2 r lV2 



1 1/2 



and 



= C||||(<7||7_l) where Cj_x, 
justable constants. Then 7^ = [C 
<r ± = [(2/3)Jd 3 kk 2 I ± /Jd 3 kk 2 ± I A 

[(1/6) J d 3 kkp\\/ J d 3 kk 2 _I 1 _] 1/2 , where we have set 
C||_L = (2/3)Cxx, C|||| = (l/6)Cxx to ensure that 
(j j_ =2/3 and cr m = 1/6 in the isotropic case. 

In order to model gradual anisotropization of the veloc- 
ity statistics by the back reaction, we define the stretch- 
ing wave number k s (t) such that the total magnetic en- 
ergy W(t) at time t is equal to the energy of the hydrody- 
namic eddies at k > k s (before they feel the nonlinearity) . 
We assume that the eddies at k < k s remain isotropic 
(unaffected by back reaction), while those at k > k s are 
two-dimensionalized. Specifically, for kg < k < k s (t), let 

47rfc 2 /«(fc, \n\) =E(k), J«(k, \fi\) = 0, (10) 



while for k s (t) < k <k v , 

47rfc 2 / (i) (fc, \n\) 
47rfc 2 / (a) (fc, \n\) 



2r 2D E{k)S(^) 7 
2E(k)5(fi). 



(11) 
(12) 



Here k s (t) is defined by c 2 £"E(k) = W(t), I® and 
/( a ) are coefficients of analogous to «W and «w 
[Eq. |JSJ], fco and k v are the forcing and viscous wave 
numbers, C2 and r 2 ^> are adjustable parameters. We take 
E(k) = C K 6 2 / 3 k~ 5 / 3 (with C K = 1.5) for fc e [*o,fe/]. 
The specific form of E(k) will only affect details of the 
transient time evolution, not the saturated state. E(k) 
will not figure in what follows, because it multiplies fiS(fi) 



in all relevant expressions. Coefficients in Eq. (JSJ) now de- 
pend on W(t): a straightforward calculation gives 



a ± (t) 

r(t) 
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4m, (i) 



(1 + VK /^) 5 



-1/2 



[rw] 1/2 ,(i3) 



(1 + Wo/W„) 2 

1 



(1+ W(t)/W,,) a 

1/2 

[r(tr 1/2 , 



(14) 



(1 + W(t)/W„) 2 (1 + W /W„) 2 



f2D 



(1 + w{t)/w„y 



where 7 = c x [/*" dk k 2 E(k)\ 1/2 , c x = [(5/18)C X x] 

the viscous-eddy energy is W v /c 2 — (3/2)CK£ 2 ^ 3 fci7 2 ^ 3 , 
and the total energy of the velocity field (before suppres- 
sion) is W a /c 2 = J k "dkE(k). Equations (|13I - I14|I rep- 
resent a generalization of the model first introduced in 
Ref. Q and reduce to it when r 2 D = 0. They include the 
effect of quasi-2D mixing of the folded magnetic fields by 
eddies whose stretching component has been suppressed. 
The spectrum of these mixing motions is modelled by 
Eq. (jll|l , where r 2 D parametrizes the strength of the mix- 
ing relative to the original unsuppressed 3D turbulence. 

The behavior of our model is easy to predict. The kine- 
matic growth stage [71 = (6/5)7, <7± — 2/3, mi = 1/6, 
and s = 3/2, 7 = (3/4)7 in Eq. ©] lasts until the total 
magnetic energy reaches the energy of the viscous-scale 
eddies, W ~ W v . After that, the velocity is gradually 
anisotropized, stretching is weakened, but mixing con- 
tinues at k > k s (t). A steady solution is reached as soon 
as g j_ and mi have decreased enough to render 7 = in 
Eq. 10. This gives a± = ±o\\ ~ 0.078. The correspond- 
ing spectral exponent in the interval k v <C k <C k v is 
s ~ 0.23. This solution is valid in the limit k^ ^ k v 
(Pr m ^> 1), but convergence in Pr m is only logarith- 
mic. In practice, numerical solution of Eq. JSJl shows 
that a scale separation of many decades is required for 
the scaling to be discernible. This is not achievable in 
direct numerical simulations. We have, therefore, solved 
Eq. JSJl with the same parameters as those used in our 
simulations There are three adjustable constants: 
ci, C2, and r 2 D- The solution does not, however, depend 
very strongly on them: c 2 is irrelevant as it amounts 
to overall rescaling of energy, r 2 D has to vary by an or- 
der of magnitude to cause significant change, and even 
ci (which affects the value of k v ) does not require very 
fine tuning. We have compared the model solutions for 
the same fixed values c\ = 1/3 and r 2 D = 4/5 with the 
(normalized) spectra obtained in numerical simulations. 
A sequence of runs with large Pr m and low Re (the so- 
called viscosity-dominated limit: a necessary compromise 
at current resolutions U) is well fitted by our model 
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FIG. 2: The bold lines are the normalized saturated energy 
spectra for simulations in the viscosity-dominated regime Q. 
The thin lines are the spectra predicted by our model. 

in both kinematic (not shown) and nonlinear (Fig. |2J 
regimes (except at k/2% = 1, 2, where finite-box effects 
are important). Note that the velocity field in these runs 
is random [because of random forcing, Eq. Q], but, un- 
like in real turbulence, spatially smooth and one-scale. 

It is extraordinary that our minimal model has repro- 
duced non-asymptotic numerical spectra so well. We do 
not claim that it constitutes a quantitative theory of 
nonlinear dynamo. It does, however, provide a simple 
demonstration that the available numerical data is con- 
sistent with magnetic-energy spectra exhibiting a very 
flat positive spectral exponent in the interval k v -c k <C 
krj if sufficiently large scale separations were resolved. 

It is clear that the viscosity-dominated simulations 
(low Re) are described very well by our model. The 
case Re > 1, Pr m ^ 1 is much harder to tackle. If 
mixing by velocities at k e [k si k„] remains efficient 
[as implied by our 2D approximation (|11H12|) ]. then k s 
stabilizes at a value ~ k v and saturated magnetic en- 
ergy scales with Re as the energy of the viscous eddies, 
(B 2 ) ~ Re _1 ^ 2 (it 2 ). This outcome does not appear to 
be borne out by the available numerical evidence, which 
rather suggests (B 2 ) < (u 2 ) (though limited reso- 
lutions preclude a definitive statement). In our runs 
with Pr m = 10 and Taylor-microscale Reynolds number 
Re > ~ 45 (Re - 100) and with Pr m = 1, Rex — 155 
(Re ~ 400), our model in its present form overestimates 
the magnetic energy at large k, but underestimates it 
at low k (Fig. ^): an indication of too much mixing 
in the model |l8|. Indeed, when Re > 1, the nature 
of the anisotropized velocities in the interval [fc s ,fc„] can 



be very different from the interchange-like motions that 
give the 2D mixing in the viscosity-dominated case. In 
Ref. 0], we argued that the interval [k a ,k u ] is popu- 
lated by Alfven waves that propagate along the folds. 
The saturated spectrum is then the result of a super- 
position of waves and folds (which accounts for the large 
amount of small-scale magnetic energy) . Since the Alfven 
waves are dissipated by viscosity, they can only exist if 
the stretching scale becomes much larger than the vis- 
cous scale: possibly as large as the outer scale (k s ~ fco, 
cf. Q). This is only allowed if the waves do not mix 
magnetic field as efficiently as the interchange motions 
do. For our model, the required modification would be 
that the mixing rate 7^ should decrease with k s . The dy- 
namo saturation would then be due to a balance between 
stretching and mixing by partially anisotropized motions 
at the stretching scale. 

Detecting Alfven waves along folds is a challenge for fu- 
ture numerical work. The main conclusion of the present 
study is that the nonlinear dynamo in a random one-scale 
flow can be described by a simple model where saturation 
is achieved via partial anisotropization of the ambient 
velocity, a result quantitatively supported by agreement 
with direct numerical simulations. 
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